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1 Introduction 



The recent experimental advances on the cold Fermi atoms [l]-[4] have revived a great 
amount of theoretical interest in superfluid fermion systems. Prom the standard Bardeen 
Cooper Schrieffer (BCS) theory of superconductivity [5] , it is known that a Fermi system 
with an arbitrarily weak attractive interaction induces a pairing mechanism, which takes 
place between the spin-up and spin-down particles with equal and opposite momenta on 
the common Fermi surface. The system, being in the fully gapped BCS state, exhibits 
s-wave superfluid properties and opens a homogeneous gap in the excitation spectrum. 
However, when the system is placed in the presence of a weak magnetic field, the Cooper 
pairing mechanism becomes unstable due to the mismatch of Fermi surfaces of the 
opposite spin particles. The pairing is now likely to occur at a nonzero total momentum 
in the so-called Larkin, Ovchinnikov, Ferrel, and Fulde (LOFF) phase, where the gap 
parameter acquires a spatial variation and breaks translational symmetry [6]. 

Fermion systems with mismatched Fermi surfaces and still having spherical symmet- 
ric interaction can be found in the models of asymmetrical two-component fermions, for 
instance, when we attempt to pair fermions with different species or attempt to pair the 
same fermion at finite polarization. A state of this type was first studied by Sarma long 
time ago [7] . The state is said in the Sarma phase. Because of the two distinct Fermi 
surfaces, the Sarma phase contains the phase separation in momentum space, including 
the paired region about the Fermi surfaces and the unpaired region in between the Fermi 
surfaces. As a result, it is frequently termed the breached pair superfluid state (BP) 
[8]. Besides, the same state is sometimes called the gapless superconducting phase, since 
a single quasiparticle dispersion gets across the free Fermi surface [9]. In general, the 
gapless state of matter is characterized by the spatially homogeneous composition of the 
paired fermions, forming a superfluid, and the unpaired fermions, retaining the property 
of a Fermi liquid. 

In the weak negative coupling regime, it is recently pointed out that the spatially 
homogeneous Sarma phase (the BP state) is unhkely to be the favored ground state for 
the asymmetrical fermion system [10]. Here, the weak-coupling regime is relevant to the 
limit of small gap parameter of the system. In the context of grand canonical ensemble, 
on one hand, it is shown [10] that either the normal or the BCS states is the energeti- 
cally preferable ground state, depending on the degree of asymmetry between the Fermi 
surfaces of the two species of fermions. On the other hand, the spatially inhomogeneous 
mixed phase consisting of BCS superfluid and normal states is suggested to be the true 
ground state in the context of canonical ensemble. Despite this, it does not rule out 
the possibility of constructing models in which the ground states are indeed charac- 
terized by a spatially homogeneous state coexisting superfluid and normal fluid. One 
model exhibiting momentum dependent interaction with a large mass ratio is proposed 
to reahze the stable Sarma phase (BP) [11]. Another model, which retains spatially 
homogeneous property but breaks rotational symmetry, is discussed in [12] using the 
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anisotropic pairing interaction in the weak p-wave coupling. 

In the strong-couphng regime, i.e., the moderate values of gap parameter, the search 
for the ground states of the asymmetrical fermion system is also reported [13, 14]. It 
is shown in the former reference [13] that a homogeneous superfluid phase with a small 
polarization (or a small asymmetry in the fermion species) and nontrivial gapless ex- 
citations may be favored as the ground state over the inhomogeneous mixed phase. 
Furthermore, at infinite coupling constant, the mixed and homogeneous phases are ac- 
tually found degenerate. That is, the gapless superfiuid can exist at very strong coupling 
and the state may be accessible in cold atom experiments. The findings are consistent 
with the results shown in the article [15]. 

The purpose of the present paper is to provide an analytical study on the asymmet- 
rical fermion system, in order to gain clear insight on the spatially homogeneous Sarma 
phase in both small and large pairing gap regimes. The system we shall consider is de- 
scribed by the Hamiltonian of two species of fermions interacting through an attractive 
contact interaction at the zero temperature. Based on the mean-field analysis, a number 
of physical quantities of the system, such as the renormalized coupling constant, the 
thermodynamic potential, and the Helmholtz free energy, are formulated into of several 
integral functions. After performing the integration, we are able to express the quanti- 
ties in terms of linear combinations of the elliptic integrals of the first and second kinds. 
Consequently, the obtained analytical results allow us to investigate the asymmetrical 
fermion system in the small and large gap parameter regimes, respectively. In the limit 
of small pairing gap, we discuss and confirm on the issues of how the asymmetry in 
fermion species affects the phases of the ground state [10]. In the limit of large pairing 
gap, we find that there are two candidates competing for the ground state. Depending 
on the degree of asymmetry, it can be either a pure Fermi fluid plus a mixed condensate 
or a pure Fermi fluid plus a pure boson condensate. Moreover, it in principle enables us 
to interpolate the system in both regimes, too. 

The plan of the paper is as follows. In Section 2 we introduce the model Hamiltonian 
and derive a number of mean- field equations. In Section 3 and 4 we solve these mean- field 
equations exactly at the zero temperature when two species of fermions have unequal 
number densities or chemical potentials. In Section 5 we analyze these equations with 
two species of fermions having the same number density. We conclude the paper in 
Section 6. 

2 The Model and Mean-Field Equations 

Let us consider a system of two species of fermions interacting through a short-range 
attractive potential, i.e., a dilute system with asymmetrical fermion species. The first 
species, denoted by A, has mass Ma and total particle number A^^. The other species, 
denoted by a, has mass Ma and total number A^a. Without loss of generality, we set 
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Ma = mMa with m > 1 and Na = nNa while keeping n arbitrary (0 < n < oo) in 
the discussions. The nonrelativistic version of the system is well described by the grand 
canonical Hamiltonian 

p p,g 
where p and q are momenta, Ap, ap, ^41 and at are the annihilation and creation op- 



erators, g is the bare coupling constant and taken to be negative, in addition = 
£A,a — fJ'A,a stand for the dispersion relations of the A-type and a-type fcrmions, re- 
spectively. Here, eA,a = /'^MA,a ai'e the kinetic energies and iJbA,a are the chemical 
potentials. The notation p = |p| is adopted. 

Following the mean-field analysis of the BCS theory, we introduce in the Hamiltonian 
(1) a real order parameter, the pairing gap, as A = —gJ2p{C'-p^p) ■ This allows us to 
simplify the Hamiltonian that can be further diagonalizcd via the Bogoliubov-Valatin 
transformation [16] . It turns out that the diagonalizcd Hamiltonian is of the form 

^ = E {K B'pBp + K hlph.p) + E {Cp -E^)-^. (2) 

p p ^ 

Here, Bp, hp, Bp, and fot arc annihilation and creation operators of the so-called B- 
type and 6-type quasi-fermions, whose corresponding energies are denoted by and 
-Ep, respectively. The quasi-fermions B and h are mixture of the original fermions A 
and a and play the roles of the quasihole and quasiparticle excitations to the ground 
state of the mean-field Hamiltonian. Explicitly, we have the operator transformation 
Bp =UpAp — Vp a^p and b-p = Up a-p + Vp A-, with the transformation functions Up 
and Vp fulfill the relations: 

ul + vl^l and ul-vl= , ^^^^^^^ (3) 

In Eq. (2), the quasiparticle energies and E'^ are found to be 
Ep ^ - 



2 

Bt = i 

2 



p 



-{^t - Cp) + ^M+W+^'] ■ (5) 



If we designate the ground state of the system by |^') [9, 10], at a finite temperature T, 
the quasi-fermions will have occupation numbers described by (\1'| B^Bp\^) = fpiE^) 
and {'^\h^_ph-p\^) — fpiEp), respectively. Here, the Fermi-Dirac distribution func- 
tion is fpiE) — -^^w^: with P — -j^. Therefore, it is straightforward to derive the 
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thermodynamic potential Q of the system at the ground state | ^) 



(6) 



The mean particle number density of each type of fermions A and a can be carried 
out by applying the thermodynamic property of the fl{A, ha, iJ^a) as 



d 



I 



d^p 



(27r)3 

d 

n{A,HA,Ha) 



d^p 



(27r)= 

Furthermore, a consistent equation for the gap parameter A is obtained by 

p 

d^p 



(7) 



(8) 



-9R J [upv, (l - MEi) - ME^)) - 



(9) 



Because of the choice of contact pairing interaction in (1), the gap equation diverges 
at the ultraviolet limits. Thus, we are lead to introduce a temperature-independent- 
regulator to regularize the gap equation in the second line of Eq. (9). The quantity Qr 
is the renormalized coupling constant that is related to the two-body (S-wave) scattering 
length in the low-density approximation by = 2{^i)-k ■ 

Besides the thermodynamic potential (6), we can as well perform the Legendre trans- 
formation to find the Helmholtz free energy T of the system 



JP(A, Na, No) = n{A, iiA, l^a) + I^aNa + A^aA^a , 



(10) 



where the chemical potentials ^a and are presumably calculated in terms of the mean 
particle numbers Na and Na by inverting Eqs. (7) and (8). 

At the zero temperature, the system is still properly described by the above equations 
(6-10), however, with the Fermi-Dirac function takes simple numbers. Depending on the 
value of energy E, we have friE) = 1 for £; < and /f(-E) = for > 0. With the 
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explicit expressions of (4) and (5), the zeros of quasiparticle energies and are 
located at the energy values 



For the purpose of simplicity, we have introduced in this equation the scaled chemical 
potential of a- type fermions as fia = ^ and the scaled gap parameter as A = 

Both quasiparticle energies E^ and E^^ are positive definite when e± defined above 
(11) are found to be complex conjugate to each other. However, when and £_ are both 
real, then one of the quasiparticle energies, either Ep or E^, can take negative values in 
the energy interval e e [£_,£+]. Note that £+ > £_ and £_ may become negative. 

We are now in a position to investigate the analytical results of the physical quanti- 
ties, such as fiA, fia, Or, ^, and JF introduced above, when the system is at the ground 
state at the zero temperature. We shall discuss three possibilities of interest below. They 
are separately the asymmetrical superfluid with Na > Na, the asymmetrical superfluid 
with Na < Na, and the symmetrical superfluid with Na — Na- 



3 Asymmetrical Superfluid with A^^ > Na 

The system is in one branch of the so-called Sarma phase (the BP phase), which can 
accommodate at the same time a Fermi liquid and a superfiuid in the weak-coupling 
regime. In this branch and thereinafter named the first branch, for a given value of gap 
parameter A, the quasiparticle energy E^ will become negative in the energy domain 
(£_ < e < £+), whereas E^ is still kept positive. From Eq. (11), it is shown that the 
allowed values of ^a and fia reside in the interval pia — ^a £ (— cxd, — 2A]. Note that 
fia = ^ and A = When Jia — ijla = —2 A, the boundary points of the energy 

interval match, that is, S- = £+. Moreover, the lower bound of the interval e_ > 
for pa £ (~^' A'-A ~ £_ < for fla G (~oO) ^^^6*^ analysis, the 

thermodynamic potential (6) can be reformulated into 



1]«(A 

mi I^Ai l^a J 

_ {2MAf'^ f 
47r2 1 




de\/e 



{e - iiT? + 2A 



^ ^ de^/e {e - /xa) | , 



^(£-//t)2 + (2A^)2 



(12) 



where the kinetic energy e = £a— ^-type fermions is used as the integration 

variable. The new notation for the gap parameter is defined by A^ = ^^A_ = ^^^A and 
the total chemical potential is jiT = ^^^^ = ° • Here, the superscript of the 



5 



Q implies that we are in the first branch of Sarma phase, such that Na > Na (see Eq. 
(14) and the paragraph below) or, equivalently, P-a^ 2 A. 

In deriving Eq. (12), we have combined the gap parameter equation (9) with the 
thermodynamic potential Eq. (6) for the elimination of the bare coupling g. Choosing 
the same integration variable £, the gap equation (9) becomes the integral equation for 
the renormalized coupling constant qr, which takes the following form 



M' 



3/2 



(m + 1) y27r2 




+ j de^/e 



In the same manner, the particle number density equations, Eqs. (7) 
reshuffled into 



NA-Na 



3/2 



2A/27r2 3 



3/2 



3/2 



1 - 



£- jJLT 



(e - fir)' + (2A^)^ 



^+ de 

(13) 

and (8), are 



(14) 
(15) 



Here, we obtain an algebraic equation plus an integral equation in Eqs. (14) and (15), 
instead of the original two integral equations given in Eqs. (7) and (8). Since e+ > e_, 
in Eq. (14) the number of A-type particles has to exceed that of a- type particles. This 
validates our earlier claim. That is to say, some of A-type particles will remain in the 
normal state even after the pairing mechanism occurs. As a result, we expect that 
the chemical potential ^ia of A-type fermions will stay in the positive region and never 
penetrate the zero value for any finite gap parameter Am- 

At zero pairing gap A^ = 0, the chemical potentials for both species of fermions 
are easily calculated. It is found that /^^(A^ = 0) = £^ and Jia^^m = 0) = (^)^'^^£^, 

where £^ = '^^^^Ma •^'^^^ Fermi energy of the 74-type fermions, that is, the 

unpaired Fermi energy. Here, n = ^ > 1 is the particle number ratio. It is pointing 
out that when the combined factor m(^)^/^ > 1 with the mass ratio m = then 
lia = infla > IJ'A at A„i = 0. As the gap parameter increases under the circumstances 
of fixing n a constant, the value of will decrease in response and eventually go to 
negative region, signifying the formation of molecular bound state. However, the value 
of /lA should remain in the positive region. We thus expect a crossing between fj,A and /la 
curves that will happen at some finite value of A^. The particular solutions of chemical 
potentials so that //^(Am) = //a(Am) can be found by solving the number equation (14). 
Generically, it is an algebraic equation of degree six and will simplify to an equation of 
degree four for the choice of mass ratio m = 1. On the contrary, when the combined 
factor m(^)^/^ < 1, we should observe no crossing between the two chemical potential 
curves /ia and /la to occur at any value of A^. 
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All the equations (12)-(15) shown above are valid only for the case when e_ > 0. 
If we notice at some gap parameter that e_ < 0, then these equations need to be 
modified by simply letting s_ — 0. It is because that negative kinetic energy is not 
physically allowed. 

To continue on the discussion of the analytical results, it is convenient to use dimen- 
sionless quantities. For this purpose, we follow and adopt closely the conventions and 
notations of the reference [17]. Let us introduce the dimensionless energy variables and 
dimensionless total chemical potential by 



2A^ 



and 



2A^ 



(16) 



and then express Eqs. (12)-(15) using them. 

In the sequel, the analytical expressions can be obtained and will be presented in 
following paragraphs. It is found that, in terms of the dimensionless quantities a± and 
Xq, the thermodynamic potential (12) becomes 



m + 1 (4M^A„)5/2 



J'\xo, 0) + a;W(a;o, a+) - u;^'\xo, a_) 



47r2 2Ma 
Similarly, the gap equation (13) reduces to 

9rMa 



[17) 



1 = 



(4MaA^)^/2 i^(^xo,0)+Ii{xo,a+)-Ii{xo,a-)-{^ - ^) , (18) 



(m + l)7r2 

and the number density equations (14) and (15) are expressible as 

1 



A^4 - iV„ 



67r2 



Na = -^{^MA^mf'lhixo, 0) + hix^, «+) - hix^, «_) 



(19) 
(20) 



Moreover, the substitution of Eqs. (17), (19), and (20) into the Helmholtz free energy 
(10) renders 



■^^^^(A^,/iA,/^a) 



m 



+ 1 (4MaA^)5/2 



47r2 



2Ma 



(^0, 0) + /W (xo, a+) - /^^^xo, a-) 



(21) 



Because the Helmholtz free energy is a function of number density, we have in mind in 
this equation that jiA = A*a(A^a, Nq) and //„ = /J'aiNA, Ng) and the relations are provided 
in Eqs. (19) and (20). 
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The functions Ii{xQ,a) and 72(2^0,0;) appearing in Eqs. (18) and (20) are integral 
functions and given by 



Ii{xo,a) = 

POO 

l2{xo,a) = / 



dxx^ 



dxx 



^{x^ - Xof + 1 



X — Xq 



x^ 



xoy + 1 



(22) 



(23) 



where dimensionless. Also, the functions u^^\xo,a) and f^^\xo,a) in Eqs. 

(17) and (21) are short-hand notation of the combinations 



1 2 

uj^^\xo, a) = hixo, oi) - -Ii{xo, a) - Xohixo, a) H — 

^ Tft n" J- 



a 



5/2 



3 2A„ 



/^^^ {xo, a) = 77(^0, a) - (xo, a) + 



2 «5/2 



where the function I^{xQ, a) above is defined as 



roo 

j{xq, a) ^ I 

J^/a 



dxx 



1 - 



X — Xq 



^(x2 - xo^ + 1 



(24) 
(25) 

(26) 



It is stressed that this integral is well defined at the ultraviolet limits. 

The detailed calculations of the integral functions Ii{xQ,a), l2{xQ,a), and Ii{xQ,a) 
are straightforward but somewhat tedious. The derivation is similar to that reported 
in the reference [17]. Therefore, we only present the results and skip the cumbersome 
parts. They can, respectively, be performed and expressed as linear combinations of the 
elliptic integrals of the first kind F(0, k) and second kind E{(I),k). To be precise, we 
obtain the following exact results 



Ii{xo,a) = \fa 



x\ 



1 



(27) 



a - J {a - xof + 11 + 



Xq 



a 



+yr+ 



^0, 



1 



Irixn, a) — 

7V 0, y 



1+xl- Xoj F{(t)a, H) + 2XoE{(j)a, n) 

5 - 2q;^ + 2^ {a - Xq)^ + 1 a + + 



(28) 



3 



a 



l + xt>, 



+ + (3 - + xo^l + xl) F{(t>^, n) + 2{xl - 3)E{(j>^, «) 



(29) 
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Here, the angle 0q, and the moduhis k appearing in the arguments of the eUiptic integrals 
F(0, k) and £^(0, k) are defined by 



(pa = COS ^ 



a 



1 + 



a 



1 + xl 



and 



K 



1 + 



Xq 



^1 + 2^0. 



(30) 



(31) 



Note that < 1. Now wc can compare the formulas of Ii{xq, a) and l2{xo, a) to those 
calculated in the paper [17]. When letting a = so that 0^ = vr, it is found that Ii{xq, 0) 
and /2(a;o,0) given in Eqs. (22) and (23) reproduce exactly the desired answers. 

Besides the analytical results, it is straightforward to check the approximate expres- 
sions of our results in the limits, when xq — > ±oo. This can be achieved by simply 
retaining a few significant terms in the corresponding expansions. 

We will analyze and present these approximate results in the case that the particle 
numbers A^^ and Na are both fixed quantities. Hence, the particle number ratio n is also 
fixed. The discussion here is relevant to the system of trapped cold fermionic atoms, 
where the atoms are confined in a real space. A remark regarding to the case when the 
total particle number is a fixed quantity is also illustrated. 

The first hmit of interest is when ^ 1 so that «;^~1 — ^ = 1 — (^) • Thus, 

it indicates the small gap approximation — > 0, or simply it means <C Sp. We 

define three constants by = 1, Oq = (^^) and Tq = . Under such definition, 

the constants satisfy the inequality < Tq < Aq. Then, via the Eqs. (19) and (20), the 
chemical potentials /x^'' and fi^^^ are found to have the following serial expansions 



/^?(A^) 



t 



m + 1 



ao 



( m + 1 




+ 



^(^o,ao) 



+ o{e) 



^ \m{^/A'Q - ^) ' 2y/To 
where t = is the expansion parameter and the ^-function is defined by 

y{AQ, ao) = In 



(32) 



(33) 



(34) 



In Eqs. (32) and (33), we observe in the limit A^ < that i^A{Am) - A^a (0) > 
and fi^K^rn) — < 0- Both the chemical potentials are adjusted accordingly to 

accommodate the change of the band structure due to the nonzero pairing gap. 
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Similarly, the gap equation (18) turns out to be 



1 



MAk 



[m + l)7r2 



Mo + TT- ^(A, oo) 



+ 



t 



8ToVao^o 



2(1 + m) (yA^ + m^)(2To- v/^^) 



VapAo - 4To 



^(>lo, ao) + 



'Aq - ^ 



m + 1 



6?(Ao,ao)^ 



+ 0(^2) , (35) 



where kp — 2MaSp — (67r^iVA)^/^ is the free Fermi momentum of the A-type fermions 



at 



0. Recall that 



Ml! 

2Ma 



is the associate free Fermi energy. 



It is apparent that the solution of Eq. (35) can exist only for ^ — ^ > 0, where we 
have suppressed the superscript in the qr. The critical coupling constant Qc is [9, 10] 



1 

9c 



MAk^ 
[m + l)7r^ 



' Ao + ■ 



^0 



Q{Aq, ao) 



(36) 



and can take positive or negative values entirely determined by the sign of the object in 
the parenthesis. In Essence, the sign of Qc depends on the complicate combination of n 
and m values, and has ^ — > — oo as the number ratio n approaches unity. 

Furthermore, the thermodynamic potential (17) has the ^-expanded series (in the 
hmit < Bp) 



+ 



Att'^Ma 
(m + 1) t 



15 




g{A^,ao)\+0{e 



(37) 



From the expression of the l^^^H^m) in Eq. (37), wc arrive at an interesting conclusion 
concerning the status of the normal state in the grand canonical ensemble. That is 
that the normal state is found stable, a global minimum, when the critical coupling 
Qc < 0, whereas it becomes unstable, a local maximum, when > 0. Numerically, if we 
choose the mass ratio m = 2, then the critical coupling reverses its sign at n ~ 4.341. 
Moreover, the Helmholtz free energy (21) is given by 



+ maf)+0(t^) 



(38) 



When setting i = 0, Eq. (38) is nothing but the Helmholtz free energy added separately 
from the unpaired A-type and unpaired a-type fermions. 
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In the same limit ^ £f, another important case to discuss is when the total 
particle number (A^^ + Na) is kept fixed, but the conversion between A-type and a-type 
particle numbers is allowed. This case is worth studying in the high energy physics. The 
overall quark number is protected by the conservation of baryonic number, but quarks 
themselves can interchange their flavors due to the weak interaction. In this case, we 
can easily extract the approximate results from Eq. (32) to Eq. (38) via the following 
simple replacements: the Fermi momentum kp is replaced by kp = {6n'^{NA + iVa))^/"^, 

/ i,T\2 

the Fermi energy ep is replaced by e'^ — , the expansion parameter t is replaced by 
t' — , the constant Aq is replaced by A'q — (;^) ^ , the constant Oq is replaced 

by a'o = (;^) ^ ) lastly the constant To is replaced by Tq = . For example, 

the chemical potential /i^'*(A^) under which (A^^ + Na) is a flxed constant has the 
^'-expansion 




(39) 



All other quantities can be obtained by the analogous procedure. 

Now supposed that we want to study how the status of the normal state evolves in 
the grand canonical ensemble, using the analytical result (17), as we vary the particle 
number difference (Na — Na) = Na{n — 1) and keep the total particle number (A^^ + Na) 
— constant. We see the following facts. At the beginning, for a small value of (n — 1), 
the normal state is just a local minimum and there exists a global minimum at some 
value of Am ^ 0. As the (n — 1) value increases, the depth of normal state grows deeper 
in the {il^^^ vs. A^) plot. Then, at some particular value of (n — 1), the system exhibits 
coexisting double minima, where both of them are stable. Numerically, it happens 
roughly around n ~ 1.23 with m = 2. At this stage, the spatially homogeneous Sarma 
phase is unstable with respect to the decay into a spatially inhomogeneous mixed phase 
[18, 10]. If (n—l) gets larger and larger, the normal state now replaces the former one and 
becomes the unique global minimum of the system. Nevertheless, for even much larger 
values of (n — 1), the normal state will turn out to be unstable eventually. It becomes 
a local maximum. Using the t'-expanded version of the thermodynamic potential il^^^ 
(37), the unstable point is found numerically at n ^ 4.431 for m — 2. 

To this end, let us discuss the other extreme of interest that is defined by the limits 

Xq <0 and \xq\ » 1, so that c::^ — (^) is a very small quantity. Equivalently, it 

stands for the situation that [sp]^ <C A^,^ <^ /if,. In such limit, £_ < always holds, we 
therefore need to set S- = in the derivation of the approximate results. In the equations 

to be presented, we denote two constants by Aq — {\ — ^ and oq = (^) > then use 
J — \ -^\ as the expansion parameter. 
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It is found that the chemical potentials in the limit — > — oo are 



fc_F 



A, 



1+ml o, 

+ —3 — - + o(r3) 



Also, the gap equation can be written by 

1 MAkf- 



«o + —7^ + 0{t- 
m 



m 



dot + 



7Ao 



3/2 



3 So 



137r 



^2 



(40) 



(41) 



(42) 



Hence, the system in the weak positive coupling regime qr — > O"*" is characterized by the 
limit of large pairing gap A„ ^ Sp. It is interesting to note the implication of the above 
equations. Since the renormalized coupling constant qr is expressible in terms of the 



two-body scattering length as gn = ^^"^^'^ o,s, we are able to verify using Eqs. (41) and 



(42) that //(i) = mfi^^ 



7.(1) 



m+1 



2(m+l)7r 
Ma 

2Mj^a1 j^^^ binding energy of the two-particle, A and 

a, bound state. In addition, from Eq. (40), the jJ^^ — the chemical potential of 

leftover A-type fermions that in fact become noninteracting after the formation of the 
tightly bound diatomic molecules. Note that the number of the leftover A-type fermions 
is {Na - Na) = Na{1-^)= NaAI^\ 

The thermodynamic potential admits the approximate form 



15 (m 



2_if 
5 Oo 



17 n 
40 So^^ 



t 



+ o{t- 



(43) 



Since the thermodynamic potential is related to the negative of the pressure of the 
system, Eq. (43) shows that the excess 74-type fermions contribute mainly to the system's 
pressure. This makes sense because all the bound-state molecules are condensate in the 
ground state, thus they barely affect the pressure of the system. 
Finally, the Helmholtz free energy is 



m 



1)(A 



A\5 
F) 



Stt^Ma 



'dot 



2A 



5/2 



+ 



/44in^^ 2971 \ 1 



5(m -|- 1) I 15 So 



20 So^^ 



t 



4) 

In the limit ^ it receives contributions primarily from the tightly bound A-a 
molecular condensate and secondarily from the leftover A-type fermions. 

Once more, if we want to express the equations (40) to (44) under the situation of 
fixed overall particle number, (Na + No) — constant, then we only need in every equation 

to replace by kj^ = {Gtt^Na + iV„))^/^ by 4, t by t', Ao by A', = {^f, 



and 



finally Oq by g'q 
Sarma phase. 



3(n + l)7r 



Let us remind that n 



^ > 1 in the first branch of the 
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4 Asymmetrical Super fluid with A^^ < A^^ 



The system is in the other branch of the Sarma phase, thus will be called the second 
branch of the Sarma phase. For a nonzero value of gap parameter A, the quasiparticle 
energy may be negative in the energy range £_<£<£+, while is still positive 
definite. Using Eq. (11), one can show that the chemical potential difference is given by 
J^a~ I^A ^ [2A, oo). The equality £_ = £+ is fulfilled when fia — I^a = 2 A. Furthermore, 
we have £_ > for jiA € l^a ~ 2A) and £_ < for jiA G (— oo, ~^)- 

In the second branch, the analytical results can be carried out in a straightforward 
manner that is similar to what we have presented in the case of the first branch, that is, 
the situation of A^^ > A^^. We find that the analytical results of the number density for 
A-type and a-type fermions, Eqs. (7) and (8), are arranged into the forms 

1 



Na-NA 
1 



x3/2 



/2(xo, 0) + /2(xo, a+) - hixo, a_) 



(45) 
(46) 



where the dimensionless quantities Xq and a± are defined in Eq. (16) and A^ = 

Since > a_, we conclude that Na > Na in the second branch. That is, some 
of the a-typc fermions will stay in the normal state after pairing occurs. The both 



chemical potentials 
and fiai^m = 0) = 



at zero gap parameter are easily derived as /^^(An 



1 A 2/3^ A 



0) 



£p. Here, n = ^ < 1. Since the combination m(^)^/^ > 1 is 



always true in the second branch, we see that > ha at = 0. As the gap parameter 
increases under the constraint that n is kept fixed, the value of /i 4 will decrease and 
eventually go to negative region, while fia {ot fia) continues to stay in positive region and 
never across the zero value. Thus, we do not expect any crossing between the chemical 
potential curves fj,A and /Xa to happen at any value of Am. 

The analytical result for the pairing gap equation in the second branch is shown to 
take the identical expression as that in the first branch. That is, the gap equation in 
this branch is nothing but Eq. (18). As for the thermodynamic potential (6), it is of 
this form 



0(2) (A 

m + 1 (4MaA, 



47r2 



2Ma 



,(2)/ 



a;'"''(xo, 0) + u^^\xo, a+) — uj^^^{xo, «_) 



(47) 



Here, the superscript ^'^^ of Vl means that we are dealing with what is defined in the second 
branch of the Sarma phase, which implies Na < Na or, equivalently, — fia ^ —2A. 
The Helmholtz free energy (10) is simplified to 



m + 1 (4MaA^ 



\5/2 



47r2 



2Ma 



/(2)(xo,0) + /(2)(xo,a+)-/(2)(^o,«_) 



(48) 
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Because the Helmholtz free energy is a function of number density, we have to solve for 
IJ'A = fJ'A^NA, No) and fia = fJ^a{NA, No) using in Eqs. (45) and (46). 

Both functions uj^'^\xo,a) and /*^^)(a;o, a) shown in (47) and (48) are of the combi- 
nations 



1 2m 
oj^'^\xo, a) = hixo, a) - -Ii{xo, a) - xohi^o, a) H — 

Zi Tfl n" J- 



3 2Ar 



(49) 



1 2m tt^/^ 

f^'^\xo,a) = Ir{xo,a) - -h(xo,a) + — — ^ • (50) 

Recall that the exact results of /-functions, Ii{xo,a), l2{xQ,a), and l7{xQ,a), are given 
in Eqs. (27), (28), and (29), respectively 

All the equations presented above are valid only when £_ > 0. For the case when 
£_ < 0, the correct equations are obtained by setting £_ = in the corresponding ones. 

Now we want to investigate the approximate results in the second branch, in the 
limits when xq ±oo. Again, we keep the number ratio n a constant in the discussion, 
that is, both particle numbers Na and Na are separately fixed. Then, we will give a 
brief remark regarding to the case of the overall number [Na + Na) fixed. 

The first interesting limit is when xq ^ 1, which signifies the small gap approximation 
— ^ or Am -C Ep, equivalently. With the help of the number equations (45) and 
(46), we find that the serial expansions of the chemical potentials in the second branch 
turn out to have exactly the same forms as those found in the first branch. Explicitly, 
we have /i^"* = ij,a^ and fi^^^ = See Eqs. (32) and (33) for the details. As far as 
the other quantities are concerned, they all take the same approximate expressions in 
both branches of Sarma phase. That is, the approximate coupling constant g^^ equals 
g]jP in Eq. (35), the approximate thermodynamic potential fi*^^^ (47) equals Q^^'^ in Eq. 
(37), and the approximate Helmholtz free energy 

(48) equals J^^^^ in Eq. (38). In 
fact, in the limit -C Sp the expansions of all the quantities in both branches match. 
However, a minor difference to be mentioned is that in the second branch the constants 
Aq, oo and Tq satisfy the inequahty ao > Tq > Aq, instead of that ao < Tq < Aq in 
the first branch. Numerically, we find in the second branch that the critical coupling 
reverses the sign at n ^ 0.019 for the choice of the mass ratio m = 2. 

With regard to the case of fixing total particle number [Na + Na), all the relevant 
equations can be deduced by the simple replacements that have been addressed in the 
context of the paragraph of the equation (39). 

That the same question arises in the second branch, similar to that in the first 
branch, is how the normal state changes in the grand canonical ensemble, when the 
number difference {Na — Na) varies and the overall number is kept fixed. Choosing 
m = 2 in the analytical expression of the thermodynamic potential (17), we find that 
the normal state will remain as a local minimum and never become a global minimum 
for the value of number ratio ranged from n = 1 to as small as n — 0.1. This implies 
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that in the second branch the system can only have a chance to exhibit coexisting double 
minima at very large asymmetry in fermion species. The specific value of n showing this 
property is roughly around 0.08. Again, if the value of n gets even small and is smaller 
than 0.019, i.e., the critical value that in the t'-expanded version reverses the sign, 
then the normal state will turn out forever to be a local maximum. 

The other interesting limit is given by Xq < and |a;o| ^ 1, which implies that 
{spY <^ <^ ji^. Again, e_ < in this limit, so that we have to set £_ = in the 

derivation. By denoting Aq — and Oq = — l) ^ , we present the approximate 

results as follows. 

It is found that the chemical potentials in the limit xq ^ —oo are 



(1 + m)io + m ao^ + 0{r^) 



(51) 



(l + m)l o' 



(52) 



(a^) expansion parameter. Also, the gap equation, which signifies 

the system in the weak positive coupling regime g^j? ^ 0"^, can be written by 



where | 



MAk^ 
[m + IW 



rj ^3/2 



IStt 



3^0 16/1, 



3/2 



(53) 



In terms of the two-body scattering length Os, we find that jji^ — ~2M^, binding 
energy of molecular bound state, and that /x^^ 



(2) 



m ao Sp 



^ is the chemical 



potential of excess a-type fermions after the formation of the tightly bound molecules. 
The excess number of a-type fermions is Na — Na = Na{ji — 1^ 

In addition, the thermodynamic potential admits the ^-expansion 



NAaT- 



A\5 



(m + l)(A;#) 



4m clq 



5/2 



2af 



ITtt 



15(m + l) I 5io t 



+ o{t- 



(54) 



Eq. (54) means that only the excess a-type fermions contribute to the system's pressure 
in the hmit | — 0, since Q = — P, the negative of the pressure. 

At last, the Helmholtz free energy in the limit ^ Sp becomes 



SnmA 



A\5 



'Aot' 



2m a, 



5/2 



44 af 



297r 



5(m + l) I 15io 20io^V ^ 



(55) 

Both the bound A-a molecular condensate and the leftover a-type fermions contribution 
to the Helmholtz free energy. 
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As regards the resembling equations when the total particle number (Na + Na) is 
fixed, the replacements below are found useful. That is, kp by kp, Sp by e'§^, t by t', do 

by a'o = and at last io by i'o = ( ^^"4^^^"' )^ A reminder is that n = ^ < 1 in 

the second branch of the Sarma phase. 



5 Symmetrical Superfluid with Na = Na 

The system is in the so-called BCS phase, in which the two species of fermions have the 
same particle number. In this phase, both quasiparticlc energies and are always 
positive for any given gap parameter A. That is, > and E^ > 0. The allowed 
chemical potentials satisfy the relation fia — l-^A ^ [— 2A,2A]. The analytical results for 
the number density of A and a fermions (7) and (8) take the simple forms 

NA^Na = ^{AMAArr^f'hixo, 0) , (56) 

where, as indicated, both A-type and a-type particles have the same number. Here, the 
number ratio n = jf- = 1. When the gap parameter vanishes, the integration of (56) is 
easily performed. It is found that /j-a = fia = for A^ = 0. Recall that A^ = 
The analytical result for the pairing gap equation (9) is written by 

1 - --^^^(4M^A^)V27i(xo, 0) . (57) 

Furthermore, the thermodynamic potential (6) and the Helmholtz free energy (10) are, 
respectively 



Q(0)rA 1 - ^+l(4M^AJ!^ 

WO). A X m + l(4M^Aj5^r 1 
•^^H^m,/^T) = —r-^ /7(xo,0) - -7i(xo, 



I-j{xo, 0) - ]^h{xQ, 0) - XQhixo, 0) 



47r2 2Ma 



0) 



(58) 
(59) 



where /Xt = ^f^i is the system's chemical potential in the BCS phase. The superscript 
means that we are investigating the quantities of the system in the BCS phase, so 

that Na = Na or, equivalently, ij,a — fia ^ [ ~2A, 2A ]. The exact results of the /i(a;o, 0), 

-^2(3^0,0), and /7(a;o,0) functions are given in Eqs. (27), (28), and (29), respectively. 
If we choose the mass ratio parameter m— 1, it is apparent that Eq. (56) to Eq. (59) 

reproduce the Leggett's equations for the BCS-BEC crossover with /i-r identified as the 

system's chemical potential and 2Am = A as the system's gap parameter [19, 20, 21]. 

For the increase of gap parameter, we expect the value of /it to decrease accordingly. 

Then, at some particular value of A, the chemical potential /it will go to negative region 

signifying the occurrence of the Bose-Einstein condensation. 
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It is easy to obtain the approximate forms in the hmits when xq 
hmit to study is for the case 3> 1, so that ^ or A 
the system's chemical potential /x^ •* has the serial expansion 



1 - U In 




±oo. The first 
■C By solving Eq. (56), 



(60) 



where t — is the expansion parameter. Similarly, the gap equation Eq. (57) is 

1 MAk 



.A 



(m + l)7r^ 



(61) 



This gives us the famous BCS value of the gap parameter at zero temperature, A^ ~ 
4e^^ e^'^/^'^-pl"*! e;^. It is clear that in the BCS phase there is always a solution to the 
gap equation, no matter how small the coupling constant gR is. This is in contrast to 
the cases discussed in the Sarma phase. 

Furthermore, in this limit A^ <^ the approximate expressions for the thermody- 
namic potential (58) and the Helmholtz free energy (59) reduce to 



(m + l)(fc# 



A\5 



4-! "7!- -'*^' 



and 



5 



(62) 



(63) 



When setting t = 0, Eq. (63) is nothing but the Helmholtz free energy of the unpaired 
A-type and unpaired a-type fermions. 

Using the approximate equations (38) and (63), it is easy to check in the context of 
canonical ensemble that the inhomogeneous mixed phase is indeed the favored ground 
state over the homogeneous Sarma phase in the limit of small pairing gap. For this 
purpose, let us consider a trapped system of two species of fermions A and a with 
separately fixed particle numbers Na and Na- As regards the Helmholtz free energies 
in the two branches of the Saram phase, they both take the same form as shown in Eq. 
(38), that is, J-'^^^A^) = JF(^)(A„J. With respect to that in the inhomogeneous mixed 
phase, it contains a Fermi liquid state of unpaired fermions plus a fully gapped BCS 
state of paired fermions. Explicitly, for the case of n = ^ > 1, the Helmholtz free 
energy in the inhomogeneous mixed phase is given by 



{k 



F^O ) 



207r2MA 



(m + 1)(A;^ 
Sti^Ma 



(64) 
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where the first term on the left-handed side is from Fermi liquid and the second term is 
from BCS phase. The surface tension effect of the interface is not taken into account in 

the analysis. Note that Aq = ( 1 — - j and = ( - In addition, for the case of 

n = ^ < 1, the Helmholtz free energy in the mixed phase is then 



t'^^) (A ) 



A -1/2x5 



+ 



(m+ l){k 



A\5 

f) 



1 1/2 



(65) 



2/3 



where Aq — 1 and olq— (^ — 1^ ' . The direct comparsion of the Helmholtz free energies 



between the Sarma and mixed phases results in {T^^^ — 

-(2) N / 5/2 



(1) 
mix/ 



(Af-af-if)>0 

for the number ratio n > 1 and (JF*^^) — T^l^) ~ m{(ij'^ — A'^'^ — Sq^^) > for the number 
ratio < n < 1. In both branches, the Sarma phase always has the higher Helmholtz 
free energy than the mixed phase does. Therefore, in the limit -C the spatially 
inhomogeneous mixed phase is found to be the energetically preferable ground state for 
the asymmetrical fermion system [10]. 

The other interesting hmit is when < and \xq\ » 1, that is, (e^)^ ^ ^ //^. 

By denoting tq = , it is found that the chemical potential is 



A 4.2 



(66) 



where 



/ e-^ \2 

( A^j expansion parameter. Also, the gap equation can be written by 

1 MAk^ 



TT ^ 

77 V^^ - 



Stt 1 



+ o{t- 



(0) 



2MAa 



In terms of the two-body scattering length a^, we obtain that fi>p 
binding energy of molecular bound state in the weak positive coupling regime. 
Furthermore, the thermodynamic potential admits the ^-expansion 



(67) 



^ is the 



(m + l)(fc^ 



A\6 



7tt 1 



(68) 



Eq. (68) indicates that the system's pressure is extremely small in the limit | 
At last, the Helmholtz free energy is 



0. 



(m + l)(A;# 
87r2MA 



A\b 



Tot' 



237r 1 



(69) 



Here, it results solely from the tightly bound A-a molecules, the standard Bose-Einstein 
(BE) condensate. 
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We end this section by posting a frequently asked question. In the large pairing 
gap limit ^ Sp, what is the preferable ground state in the context of canonical 
ensemble for the trapped cold asymmetrical fermion system? From the phase structures, 
we know that there are two competing candidates. The first one is the Sarma phase 
consisting of a pure Fermi fluid and a mixed condensate, which is characterized by the 
homogeneous mixture of fermions in a boson condensate. The other one is the separated 
phase composed of a pure Fermi fluid plus a standard BE condensate. We resolve this 
problem by computing and comparing the Helmholtz free energies of the system at the 
two phases, based on the analytical results presented in this and previous two sections. 

First, we consider a trapped system of A- type and a- type fermions with particle 
numbers A^^ > iV^. The numbers and Na arc both flxed quantities. The Helmholtz 
free energy J^(i)(A„) in the Sarma phase is given by Eq. (44). In the separated phase, 
the system is made of a standard BE condensate of tightly bound A-a molecules 
and a Fermi fluid of {Na — No) excess A- type fermions. With the help of Eq. (69), the 
Helmholtz free energy the separated phase is therefore 



•^sepl^m) 



l/2^5 

) 



TT ^ 2 237r 1 
2 20 So^^ 



- + c»(t-^) 



(70) 



where the first term on the left-handed side is from Fermi liquid and the second term 
is from BE condensate. To construct (70), we neglect the surface tension effect of the 
interface between the two components. Then, the subtraction of the Helmholtz free 
energies renders 



^i^(A„) 



Sep v 



AAA, 



3/2 



15 qq 



37r 



10 a: 



-3/2 



(71) 



Hence, we have (JF(^) ^ ^^^p^ ^ ^ ^..^ ..^...^^.x .^^.^ ,^ ^ 

1 < n < ||. Recall that = (^f^) • This implies that at large degree of asymmetry the 
system favors the Sarma phase as the ground state, while at small degree of asymmetry 
it prefers in the separated phase, the state containing a pure fermion region and a pure 
boson region. This makes sense since the more amounts the excess fermions are in the 
trap, the better chance they have to be in the boson condensate due to the degeneracy 
pressure among them. In addition, the Sarma phase and the separated phase are in fact 
degenerate in the limit t ^ oo. 

Next, we discuss another trapped system of 74-type and a- type fermions with particle 
numbers Na > Na- In the Sarma phase, the Helmholtz free energy ^(^)(A„) can be 
found in Eq. (55). While, in the separated phase, the system is made of a BE condensate 
of Na tightly bound A-a molecules and a Fermi fiuid of (Nu — Na) excess a-typc fermions. 
Based on the Eq. (69) in the BCS phase, the Helmholtz free energy J^^^^ in the separated 



scp J 



< for the number ratio n > and (JF^^) 



scp) 



> for 
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phase is now written by 



[m 



237r 1 



+ 0(0 



Similarly, we omit the surface tension effect on the interface. Note that 
Then, the subtraction of the ^^^^ and J^^^^ shows that 



44 gf 
15io 



lOio^' 



(72) 
(73) 



We have (J^^^^ -J^ie?) < for the number ratio < n < f , and (JP'^^) -T^^l) > for § < 



Sep; 



n < 1. The same previous finding applies to this second case. In the large gap parameter 
limit, the Sarma phase supports the system at large degree of asymmetry, whereas the 
separated phase supports it at smaller asymmetry. Both phases are degenerate at infinite 
pairing gap parameter. 



6 Concluding Remarks 

Based on the mean-field theory, we calculate and present the analytical results for the 
asymmetrical fermion system with attractive interaction at the zero temperature. Gener- 
ically, the results are found exact and expressible as linear combinations of the elliptic 
integrals of the first and second kinds. We discuss in details the approximate expressions 
of the analytical results in two limits: one is at small gap parameter <^ Sp and the 
other is at large gap parameter A^ » Sp. For the case of small pairing gap, we discuss 
how the normal state of free fermions evolves, depending on the asymmetry of fermion 
species. Our results also indicate that the spatially homogeneous Sarma phase is indeed 
unstable with respect to the decay into the spatially inhomogeneous mixed phase. In 
the limit of large pairing gap, we show that two candidate phases are competing for the 
ground state of the system. The Sarma phase containing a pure Fermi fluid and a mixed 
condensate is favored at large degree of asymmetry. The separated phase consisting of 
a pure Fermi fluid and a boson condensate supports the system at smaller degree of 
asymmetry. 

In addition, from our calculation it is easy to demonstrate that the existence of the 
spatially inhomogeneous mixed phase is always possible. Recall that the thermodynamic 
potential Q is a function of the variables fia, and A^ and take different forms in 
different regions of chemical potentials. In fact, the thermodynamic potential is of the 
combined expressions, Q = ^2(1) (17) in the region (j^a — fJ^a > 2A, Q = Q^^'' (47) in 
the region ha — < — 2A, and Q, — Q^^^ (58) in the region — 2A < ha — fia < 2A, 
respectively. Therefore, if we plot the (Q vs. /Za)-diagram for various values of /ia while 



20 



keeping the gap parameter fixed, we observe the foUowings. The profile of f2 is a 
continuous convcave-down function of the chemical potential /2a, and is not differentiable 
precisely at the two points, fia = fJ'A — 2A and /i„ = /i4 + 2A. As a result, the cone of 
tangents at each point describes various possible mixed phases that are composed of the 
inhomogeneous mixture of the normal Fermi liquid and the fully gapped BCS state [11]. 

Despite we have provided in the paper the analytical results for the system of asym- 
metrical fermion species with attractive interaction, we intend not to address the issues 
relating to the system in the strong-coupling regime. It is because that the Hamilto- 
nian (1) requires to obtain the analytical results might be too simplified to describe the 
complicated physical systems. To learn the system in the strong-coupling regime, one 
needs the more realistic Hamiltonian, such as at nonzero temperature, with finite-range 
interaction, and especially by the nonperturbative method, instead. 

We conclude the paper by noting that we have studied the analytical results for 
a number of physical quantities defined at the mean-field level, but we have not yet 
evaluated other physical quantities that are defined beyond the mean-field theory. For 
example, the introduction of Gaussian fluctuations on the mean-fleld background enables 
us to study the spatial fluctuations of the pairing gap parameter and the Goldstone mode 
of the broken symmetry [17, 21]. We plan to investigate the analytical results for these 
physical quantities of the asymmetrical system in the future. 

The author is grateful to D. J. Han, C. R. Lee, F. K. Men, C. H. Pao, Y. C. Tsai, 
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by Taiwan's National Science CouncU Grant No. 93-2112-M-194-010. 
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